Mesoscopics and fluctuations in networks 
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We describe fluctuations in finite-size networks with a complex distribution of connections, P(k). 
We show that the spectrum of fluctuations of the number of vertices with a given degree is Poissonian. 
These mesoscopic fluctuations are strong in the large-degree region, where P(k) < 1/N (N is 
the total number of vertices in a network), and are important in networks with fat-tailed degree 
distributions. 
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Fluctuations in finite systems, which vanish in the in- 
finite system limit, is a basic topic of mesoscopic physics. 
The study of these fluctuations is an exciting field of con- 
densed matter physics where small objects attract 
much attention. 

Despite the fast recent progress in the statistical 
physics of networks O^M0,R|, the finite-size effects in 
networks are surprisingly poorly studied. However, real 
networks, as a rule, are small objects, and their finitcncss 
is a factor of primary importance. Impressive results ob- 
tained for infinite nets often turn out to be incorrect if 
we take into account the finiteness of networks (see de- 
tailed discussion in Refs. HI?)). F° r example, the abso- 
lute random-damage stability of infinite networks with 
fat-tailed degree distributions and the absence of an epi- 
demic threshold in these nets (that is, their absolute dis- 
ease vulnerability) vanish in finite-size networks 

The finiteness of networks with fat-tailed degree distri- 
butions cutoffs the tails of the degree distributions (see, 
e.g., Refs. [p|, ^0|Jll| l ) . In this communication we consider 
another mesoscopic effect which is clearly observable in 
empirical data. We describe fluctuations of the number 
of vertices N(k) with a given number of connections k 
in a network, k is called degree, or, not quite rigorously, 
connectivity of a vertex. We consider a simple equilib- 
rium situation, where the total numbers of vertices N 
and edges L in a network are constant, and comment a 
more complex case of fluctuations in growing networks 
(see calculations in Ref. [jTo|). 

Our general results are formulated as follows: for any 
k, the fluctuations of N(k) are described by the Pois- 
son distribution. These fluctuations are of primary im- 
portance in the large-degree range, where the average 
(N(k))<l. 

A random network is, actually, a statistical ensemble G 
which includes numerous members g g G 11 1^| . These 
are different particular realizations of the network. Each 
member g enters into the ensemble with its statistical 
weight 11(g). 

When empirical researchers study a single realization 
of an ensemble, they measure N(k), i.e. the number of 
vertices of a degree k in this realization. If they have 
a possibility to investigate the complete ensemble, they 



can obtain the distribution of the fluctuations of N(k) 
at any k, that is V(N(k)). In particular, the first mo- 
ment of this distribution provides the degree distribution: 
(N(k)) = NP(k). So, studying only one member of the 
ensemble or a few ones, one can find the degree distribu- 
tion only approximately. 

There is a specific network construction, where fluc- 
tuations of N(k) are absent by definition. These are so 
called "labeled random graphs with a given degree se- 
quence" |l3| from graph theory, which are, loosely speak- 
ing, the maximally random graphs with a given sequence 
{N(k)}, ^2 k N(k) — N. Consequently, in this situation, 
P(k) — N(k)/N. These graphs are extensively used in 
numerous applications [ p^|JT^p^ , p7[ | , and, of course, the 
absence of fluctuations is related only to their N(k). So 
that, this is an exception. 

Here we study fluctuations in the dynamically con- 
structed networks that are equivalent to the random 
graphs with a given degree sequence in the thermody- 
namic limit, where fluctuations vanish. Let us consider 
an equilibrium random network with a fixed number of 
vertices N in which all vertices are statistically indepen- 
dent. Let a full set of probabilities {P(k)}, k = 0, 1, 2, . . . 
that a vertex has degree k be given, but the total number 
of edges in particular members of this ensemble be not 
fixed. Taking into account the statistical independence of 
vertices we arrive at a standard combinatorial problem. 
So, a probability that a sequence {N(k)} (J2k N(k) = N) 
is realized in a member of the ensemble has the following 
form: 



V({N(k)}) = NlH 



[P(k)] 



AT(fc) 



N(k)\ 



(1) 



This is a polynomial distribution from which the bino- 
mial form of the distribution of the numbers N(k) (the 
fluctuation spectrum of N(k)) directly follows: 

V(N(k)) = [ N N {k) )p(k) N(k) [l ~ P(k)f- N ^ , (2) 

This form, for sufficiently large networks, results in the 
Poisson spectrum of the fluctuations: 
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V{N(k)) = e -W fe )> 



(N(k)) 



N(k) 



N(k)\ 



which in turn approaches the Gaussian form 



V{N{k)) 



1 



exp 



[(N(k) - (N(k))} 
2(N(k)) 



(3) 



(4) 



when (N(k)) > 1 (or P(k) > 1/N) while the ratio 
[N(k) - (N(k))]/ v / {N(k)) is fixed, that is when relative 
fluctuations are small. Note that these general results 
are valid both for directed and undirected networks. 

Now, let the number L of edges in the realizations of 
an ensemble be fixed. In this event, the combinatorics at 
finite iV turns out to be slightly more complex, and we 
use a more convenient approach for sufficiently large net- 
works. We consider fluctuations in the dynamically con- 
structed equilibrium networks (with the fixed number of 
vertices and edges) that are equivalent to the above net- 
works with a fixed L (and to the random graphs with 
a given degree sequence) in the thermodynamic limit, 
where fluctuations vanish. Let us remind the construc- 
tion and a standard statistical mechanics formalism fl2j l 
which allows easy calculation of averages. 

(1) The ensemble includes all graphs with a given num- 
bers of vertices, N, and edges, L. 

(2) The stationary statistical weights of the realizations 
are a limiting result of the following process. Ran- 
dom ends of randomly chosen edges are rewired to 
preferentially chosen vertices at a rate /(h), where 
the rate depends on the degree ki of a target vertex. 

The function f(k) and the mean degree k = 2L/N 
determine the structure of this random network. (For 
brevity, here we consider undirected networks and do not 
discuss multiplication of edges. This will not influence 
our results.) The resulting statistical weights in this en- 
semble are products H(g) <x Y\f = iP{ki), where 



fc-i 

p(k > 0) = J] f(q) , p(0) = l 

9=0 



(5) 



At large N, the asymptotic expression for the partition 
function Z — J2 g eG ^(ff) of the ensemble is 



Z({p(k)},N,L) 



N 



where 



(6) 



(7) 



and x s is given by the equation k = x g <f>'(x s )/<S>(x s ). 

With this partition function, the following standard 
statistical mechanics relations are valid: 



d(fc) = (N(k)) = 



6\nZ({p(k)},N,L) 



5 lnp(fc) 

C 2 (k,k') = (N(k)N(k')) - (N(k))(N(k')) 

S 2 In Z({p(~k)},N,L) 
S lnp(fc)(5 lnp(fc') ' 



(8) 



and so on for the higher cumulants C n {k\, . . . k n ) which 
contain higher order products of N(k). 

Substituting Eqs. (|) and (g) into Eq. (|) readily 
shows that at sufficiently large N, when (N(k)) <C N, 

(1) C n {k\, . . . k n ) is nonzero only if k± = &2 = • • • = k n , 
so there is no correlation between numbers of vertices 
of different degrees, and that 

(2) all the cumulants are equal: C„>i(fc) = C\(k) — 
(N(k)) = NP{k). 

The latter is the characteristic property of the Poisson 
distribution. Consequently, the spectrum of the fluctua- 
tions of N(k) is Poissonian, and we again arrive at the 
same expressions @ and (|J) as for networks with fluc- 
tuating L. This universal result does not depend on the 
specific parameters of the model which can be easily gen- 
eralized. 
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FIG. 1. Typical (schematic) log- log plot of the num- 
ber N(k) of vertices of degree k in a scale-free network 
(one realization). 7 is the exponent of the degree dis- 
tribution. The fluctuations are observable in the region 
(I/7) In AT <lnfc< [1/(7- l)]lniV. 

The fluctuations of N (k) vanish in the thermodynamic 
limit if k is fixed. On the other hand, the fluctuations are 
strong in the region of large degrees where (N(k)) < 1, 
that is P(k) < 1/N. The condition P(k m ) = 1/N deter- 
mines the lower boundary k m of this region. The position 
of the cut-off of the fat-tailed degree distribution may de- 
pend on many factors. In general terms, the maximum 
possible degree of the cut-off k c is determined by the finite 
size of the network: N J, dkP ao (k) ~ 1, where Poo(k) 
is the degree distribution of the corresponding infinite 
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network (note the discussion of the complex problem of 
the cut-off position in Ref. jy}). k c is an estimate of the 
maximum vertex degree in a network and so is the upper 
boundary of the region, where the strong fluctuations 
are observable. In so-called scale-free networks, where 
Poo(k) oc fc~ 7 , the strong fluctuations should be seen in 
the broad region (l/j)lnN < lnfc < [1/(7-1)] In AT (see 
Fig. 0). 

One should note that in equilibrium networks, the fat- 
tailed regime can be realized only starting from some 
threshold value of the mean degree |h|. Below this 
threshold, a size-independent cut-off of the degree distri- 
bution sharply narrows the region of strong fluctuations. 

The Poisson spectrum of fluctuations of N(k) in equi- 
librium networks is a natural consequence of the statisti- 
cal independence of their vertices. In models of growing 
networks (e.g., see Refs. Mjl9J2Q]), the fluctuations of 
N(k) were calculated only in a small degree region and 
were found to be Gaussian [flo|| . The interesting region 
of strong fluctuations in growing nets (and in the Simon 
model is not studied yet. Although the problem is 
open, we suggest that the complete fluctuation spectrum 
of growing networks does not differ essentially from the 
spectrum of the fluctuations in equilibrium networks. 
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Ir summary, we have shown that the fluctuation spec- 



trum of the number of vertices of a given degree is 
Poissonian in equilibrium networks. These fluctuations 
are a prominent mesoscopic effect and are seen in the 
large-degree region of empirical degree distributions. As 
a rule, empirical researchers try to avoid fluctuations by 
passing to cumulative distributions. But fluctuations are 
interesting, and the simple fluctuations of N(k) that we 
have considered is only a very particular type of fluctu- 
ations in networks. We hope to attract the attention of 
empirical researches to this problem. 
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